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Closed Eigenvectors and/or Eigenvalues by Jacob i 

Method (SADOI or DOl) 

Closed Routine 

185 

thru 12 and 2n or kn locations given by pre-set 

parameter S3. -2n are required to obtain eigenvalues 

only; kn if eigenvectors are also desired, n is the 

order of the matrix. 

A block of 65 n words on the drum, given by pre-set 

parameter Sk, is used to store the given matrix. 

If the eigenvectors of the given matrix are desired, the 

block of 65(n-l) + n words immediately following this 

block is also used. 

Approximately 23n msec, for eigenvalues, 46n 5 msec. 

for both. 

n (=order of matrix) < 65. 

S2, S3, S4 



2 OOF OO cF 

3 OOF OO bF 
k OOF 00 aF 



c ia location of the Square Root 
Routine Rl. 

b is the location of 2n or kn temporary 
locations in the Williams Memory. 

a is the location of the matrix on the 
drum. 



DESCRIPTION: 



Let A be a real symmetric matrix of order n < 65 stored 
on the drum at a as described under Matrix Storage below 
(a is given by pre-set parameter 34). The sum of the 
squares of the elements of A must be less than l/2 . Let 
the Square Root Routine Rl be located as given by pre-set 
parameter S2. Then this toutine located at q when entered by 



p 


xO 


nF 




50 


PF 


p+1 


26 


qF 




— 


— 



Kite, compute the eigenvalues of A if x=5 and will compute 
the eigenvalues and eigenvectors if x=J. The mathematical 



•.« • 



method used is described in the write-up of routine 
MO. The eigenvalues H^ will be placed in the storage 
locations originally occupied by the diagonal elements 
of A, namely at a + 66k, h-O, ... , n-1. The eigen- 
vectors of A, if computed, will be stored as column 
vectors at locations beginning at a ♦ 65n. That is, 
the elements of the eigenvector corresponding to X 
will be stored at locations * ♦ 65n *- k, a ♦ 6£(n+l) +k, 
... , a ♦ 6$(2n-l) +k for k*0, ... , n-1. 
If it is desifced to obtain, instead of the matrix E of 
eigenvectors, a matrix BE, then 
(1) word 122 of this routine should be changed to 
22 103L OOF, 

<2) before this routine is entered, the elements of the 
matrix B should be stored in the locations corresponding 
to the respective elements of E (the element B.. of 
the matrix B is stored at location a.+ 6£(n+i) *j 
for i, j - 0, ... , n-1), and 

(3) this routine should be entered with the order pair 
JO nF & pF in the Q register. 

(See second and third paragraph under Matrix Storage 
below). 
MATRIX STORAGE* In order that reading from and recording on the drum 

always be done with minimum access time the first 
element of each row is stored at an interval of 65 
locations from the first element of the preceding 
row. Since the matrix A is symmetric, only the 
elements on or below the main diagonal need be stored. 
Hence, the element A ± . of the matrix A is stored at 
location a ♦ 6£L +j fpr i-0, ... , n-1 and < j £ I. 
This routine contains within it two subroirtines which 
transfer columns of elements to the memory from the 
drum and vice versa . Since these might be helpful 
to the user if used directly they are described in 
some detail in the paragraph below. 



-3- 

Iat b be the address given by pre-set parameter S3. 
I*t c be the right hand address of word Ultf, 
(relative to this routine). Let a be the drum 
address given by pre-set parameter SU. If word 
181L (relative to this routine) is the number ax 2" 39 
then the subroutine beginning at 132L (relative to 
this fcbutine) when entered using 

xO kF 

50 pT 



P*l 



26 132(M23) 



locations a ♦ 65k, a + 65k +1, ... , a + 66k -1, 

a ♦ 66k, * ♦ 66k + 65, a ♦ 66k + 130, ... , 

a * 66k ♦ 65(m»k~l) in that order to m consecutive 

locations in the Williams memory beginning at b or c 

depending qpon whether X ■ J or X - 5, respectively. 

Under the seat conditions when entered in a similar 

manner the subroutine beginning at 158L (relative 

to this reutine) dees the opposite; that is, it 

transfers words from locations beginning at b or c 

depending upon whether X - J or 5 to the drwi 

locations indicated above. Thus for k»0, 1, ... , n-1 

in that order, if the elements a. . (o < i < k) of the 
vth ki - - 

k row of the matrix A are stored in locations 

b+i and if location 181L contains n x 2~ 39 , then the 

subroutine at 158L when entered with JO kF 50 pF 

in Q will store those elements correctly on the 

drum. After routine M23 has been entered m-n or 2n 

depending upon whether only eigenvalues or both 

eigenvalues and eigenvectors were computed, and c ■ b+m. 

Hence, entering the subroutine at 132L as indicated 

±.y. 

above will cause the k eigenvalue to be stored at 
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b ♦ k or c ♦ k and (if computed) the elements of the 

th 

k eigenvector to be stored at n locations beginning 
at b f nor e ♦ n, respectively. The subroutine at 
132L uses only location as temporary storage; the 
subroutine at 1$8L uses and 1. 
ARITHMETIC CHECK* The arithmetic check described in the write-up of 

routine MO is used by this routine. If this test 
is failed the difference between the new and old 
value of the sum of the squares of the elements is 
punched out in sexadecimal form. Then an FF stop 
from the right hand side of word 112 will occur. 
This failure can be caused by a drum, memory, or 
arithmetic error. The IT stop may be by-passed with 
the white switch. 
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00 K(M23) 









KS f 








1*2 108L 




Set return address 


1 


ho hF 
U6 1*L 






2 


LU S3L 
U6 116L 






3 


10 20F 
1*2 117L 






h 


00 3F 
SS F 






$ 


1*2 183L 
10 3F 






6 


36 7L 




5 or J? 




00 IF 




J f double 


7 


U2 181L 
LU 117L 






8 


U2 117L 
U2 127L 






9 


15 37L 
Ui 181L 






10 


1*2 11UL 
1*2 131L 






11 


00 2CF 
1*6 113L 






12 


la UF 






13 


26 91L 




Iteration 


la 8F 




Ul 9F 






1U 


27 87L 
1*7 80L 






15 


& F 
$0 1% 




address ■ n 


26 


26 132L 
U7 17L 




To read drum subroutine 
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17 


JO F 
50 17L 








18 


26 132L 
L5 f 




To read drum subroutine 




19 


LU 37L 
!j2 23L 








20 


hZ 79L 

Ik 181L 








21 


ij2 22L 
U2 76L 








22 


l<9 12F 
L5 f 








23 


10 UF 
IS F 




v 




2U 


IiO IF 
80 IF 








2$ 


1<0 3F 
L3 IF 








26 


32 801 
$0 IF 








27 


L5 9F 
7U1F 








28 


LU 8F 
llO 8F 








29 


S5 F . 
ilO 9F 








30 


16 F 
bO 5F 








31 


LO I4F 
I4O 2F 








32 


ia F 

50 F 








33 


L7 3F 
L2 2F 
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M23 



3k 
35 
36 

37 
38 

39 

bo 

U* 
U5 
1*6 

Iff 
UB 

U9 
50 



36 36L 
LI IF 
66 2F 
2L 38L 
LI 2F 
10 IF 
66 3F 

85 S3 

1*0 7F 
S5 F 
hO 6F 
7J 6F 
LU U9L 
50 liOL 
26 S2 
hO 10F 
L7 7F 
10 IF 
66 10F 
KJ F 
36 U7L 
L5 115L 
£jO IF 
ljO 2F 
L5 6F 
22 57L 
SJ F 
50 U7L 
26 S2 
L5 7F 
LO 12F 
36 5UL 
L5 7F 
36 52L 



To Rl, Square Root 



To Rl, Square Root 



LOCATION 


. ' 1/ 

OH)ER 


NOTES PAGE h 


5l 


LI 119L 
22 52L 






52 


19 IF 
66 10F 






53 


S5 S3 








26 56L 






5U 


5b 6F 
75 IF 






.** 


00 IF 
LO 6F 






$6 


66 2F 






57 


S5 S3 
I4O 7F 






58 1 


L5 llUL 
22 67L 






59 


50 2F 




5 




V F 




. 


60 


J# 10F 
50 7F 






61 


7J F 








50 F 




• 


62 


Lk 10F 
IjO F 






°3 


7J 7F 
llO 10F 






& 


50 F 
36 65L 




' 


65 


7J 2F 








LO 10F 




. 


66 ; 


1*0 F 
L5 61L 






67 


fU 118L 








IrO 61L j 
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68 


U2 59L 
U2 62L 






69 


k6 6UL 
U6 66L 






70 


LO 113L 
36 59L 






71 


50 UF 
7 J IF 






72 


itO 10F 

SO IF 






73 


79 # 
LU SF 






7U 


LU 10F 
UO 10F 






75 


50 3F 
7J 6F 






76 


LU 10F 








UO F 




k * c ♦ S3 


77 


JjO 10F 

L5 Uf 




■ 


78 


LU 5F 
LO 10F 






79 


UO F 








Ul F 




k ♦ S3 


80 


J2 F 
50 80L 






81 


26 158L 
L5 17L 




To wite dm* subroutine 


82 


LU 118L 
U6 80L 






83 


U6 17L 
LO 15L 






8U 


32 17L 
36 85L 








1 
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85 


52 F 
50 85L 






86 


26 158L 
L5 15L 




To write drum subroutine 


87 


LU 118L 
1|6 15L 






88 


I46 85L 

LU 57L 






89 


1*5 30L 
U6 79L 






90 


I/) 116L 
32 1UL 






91 


Ul 7F 
U6 95L 






92 


l5 iiUl 

22 93L 






93 


F5 6F 
!0 6F 






9!* 


U2 96L 
1|2 97L 






95 


JO F 
50 95L 






96 


26 132L 
50F 




To drun reader 


97 


L5 12F 
7UF 






98 


LU7F 
UO 7F 






99 


S5 F 

Uo 5F 






100 


F5 97L 
liZ 96L 






101 


U2 97L 






l — 


LO 117L 
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102 


32 96L 
LI 11F 










103 


32 121L 
L5 95L 










lOli 


Lk 118L 
U6 95L 










105 


LO 15L 
36 93L 










106 


L5 7F 
50 11F 










107 


1*0 11F 
KD F 










108 


36 120L 










109 


22 F 




exit M23 






LO 8F 








1*0 F 










110 


V> 183L 
L2 F 










111 


36 13L 
L5 F 




• 






112 


82 UGF 
FF F 










113 


LJ F 
26 13L 




2n ♦ 3(S3) or n + 2(S3) 






21k 


V S3 










115 


50 F 

7L ii095F 

LL 1j095F 




c-2n+ 3(S3) or n * 2(33) 






116 


85 F 
00 F 










117 


75 12F 








- 


118 


7U F 

00 IF 
00 F 




(2 or 3) x (n * S3) 
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119 


20 F 
00 F 






120 


SI F 
36 13L 






121 


26 18UL 
L5 UF 






122 


32 103L 
IS 117L 






123 


U2 12UL 
L5 95L 






12U 


U6 129L 

la F 






125 


F5 12IiL 
\i2 12UL 






126 


LO 131L 
LO 181L 






127 


32 12UL 








UO F 




(2 or 3)(n + S3) 


128 


F5 127L 
UO 127L 






129 


J2 F 
SO 129L 






130 


26 158L 
22 103L 




To write drtm subroutine 


131 


N6 129L 






132 


ia F 




2n +■ 3(S3) or n ♦ 2(S3) 
Read 6nm subroutine 


K5 F 




1*2 15UL 




Set return address 


133 


32 13l*L 




5 or J? 




L5 llU 




J 


13U 


26 135L 
L5 37L 






135 


la 1U5L 








m i8il 
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136 


>2 156L 

00 9F 










137 


01 10F 
1*0 F 


1 


Save k 






138 


JO 155L | 










00 6F 




x 6b 






139 


Lb F 
Lb 157L 




*k - 65k 

•k drum read order 






UjO 


bO ll*l*L 
Ll* F 




Set to read matrix 






ii*i 


llO li*9L 
L5 ll*5L 




66k 






11*2 


Ll*F 
& 156L 










1U3 


Ll F 
36 ll*9L 










ll*b 


85 11F 
OOF 




^ead drum 






11*5 


50 F 
ItO F 




b or c 

• 






11*6 


F5 ll*i*L 
1*0 ll*l*L 










11*7 


F5 1U5L 
1*2 ll*5L 


. 








11*8 


LO 150L 
36 IbbL 










11*9 


85 11F 
00 F 




66k 






150 


JO F 
1*0 F 




(b or c) + k 






151 


F5 ll*9L 
Lb 155L 










152 


hO lb9L 








1 




F5 150L 


1 
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LOCATION 
153 
15U 
155 
156 
157 

158 
159 
160 
161 
162 
163 
16U 

165 
166 
16? 
168 
169 



ORDER 

i|2 150L 
10 156L 
36 1U9L 
22 F 
00 F 
00 6bF 
50 F 
ho F 
85 11F 
00 Sk 



NOTES 



PAGE IX) 



*3 



K5 F 
h2 180L 
32 160L 

IS ilia, 
26 161L 
L5 37L 
l* 170L 
L5 170L 
Lb 181L 
liO IF 

00 9F 

01 10F 
llO F 
JO 181L 
00 6F 

lUf 

Lii 182L 
hO 171L 
LU F 

to i76i 

L5 170L 
Ih F 
1*2 175L 
U F 



Leave drum reader 



c+c or c+b 

Matrix read constant 

Write drum subroutine 
J or 5t 



k 

x 6U 

♦k * 65k 

+ SU 

Set matrix write 

♦ k 

-66k 
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170 


32 17ft 
L*F 






171 


86 11F 
00 F 






172 


FS 171L 
li0 171L 






173 


F£ 170L 
U2 170L 






17U 


10 17ft 
32 170L 






175 


S2 17ft 
l£F 






176 


86 HF 








00 F 




66k 


177 


F$ 176L 
LU l*ft 






178 


kO 176L 
F5 17ft 






179 


U2 17ft 
LO IF 






180 


32 17ft 








22 F 




Le«r« droi writer 


181 


00 F 








00 F 




c - 2(iwS3) or (n ♦ S3) 


182 


86 UF 
00 Si* 






183 


00 I 






18U 


00 F 




8(n + S3) 


m 7F 




26 109L 






1 




1 
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